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Joseph  Kaye,  Joseph  H.  Keenan,  George  A.  Brown,  and  Robert  II.  Shoulberg. 

SUMMARY 

Reliable  experimental  data,  obtained  at  relatively  low  cost,  are 
presented  in  the  form  of  heat-transfer  coefficients  for  air  moving  at 
supersonic  speeds  in  a round  tube.  These  data  are  analyzed,  interpreted, 
and  compared  with  available  data  in  the  literature. 
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The  experimental  local  heat-transfer  coefficients  are  for  laminar, 
transitional,  and  turbulent  boundary  layers.  The  data  for  a laminar 
boundary  layer  are  given  in  this  paper,  and  the  remaining  data  will  be 
given  in  a separate  paper.  The  experimental  data  for  17  runs  are  given 
here  for  Mach  numbers  at  tube  inlet  of  2.8  and  3.0.  The  range  of  values 
of  diameter  Reynolds  number  covered  is  from  20,000  to  100,000  for  these 
laminar  flow  tests,  while  the  length  Reynolds  number  extends  to  about 
4,000,000.  The  computed  quantities  are  obtained  on  the  basis  of  a simple 
one -dimensional  flow  model,  but  a subsequent  paper  will  analyze  the  same 
data  in  greater  detail  on  the  basis  of  a two-dimensional  flow  model. 
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NOMENCLATURE 


A - cross-sectional  area,  ttD^/ 4 

A1  - heat-transfer  area.  jtDAL 

Cp  - specific  heat  at  constant  pressure 

cv  - specific  heat  at  constant  volume 

cw  - discharge  coefficient  of  nozzle 
D - inside  diameter  of  pipe 

g - acceleration  given  to  unit  mass  by  unit  force 
G - flow  per  unit  area,  w/A 

h - coefficient  of  heat  transfer,  q/A'  (t  - t ) 
k - ratio  of  specific  heats,  Cp/cv 
L - distance  from  end  of  curved  contour  of  nozzle 

M - Mach  number,  V/  v'fjkR’f 

n - summation  index,  Eq.  (8) 

Nup)  - diameter  Nusselt  number,  hD/Xm 
Nul  ~ length  Nusselt  number,  hL/Xm 
p - static  pressure 

a - rate  of  heat  transfer 

r - recovery  factor,  (taw  - tm)/(toi  - La) 

R - perfect -gas  constant 

Rejy  - diameter  Reynolds  number,  DG/pg 

Re^,  - length  Reynolds  number,  LG/pg 

S'l  - Stanton  number,  h/opG 

t - temperature,  deg  F 

T - temperature,  deg  F ubs 

V - velocity 

w - mass  rate  of  flow 

p - density 

p - viscosity 

X - thermal  conductivity 

Superscript  * refers  to  throat  of  supersonic  nozzle  w’ee;  ■;  M = 1 
Subscripts: 

aw  - adiabatic  wall  conditions 

j - station  numbers 

m - mean  stream  conditions 

0 - hypothetical  entrance  plane  of  the  tube,  where  the  boundary 

layer  is  of  zero  thickness 
of  - downstream  stagnation  conditions 

01  - upstream  stagnation  conditions 

oj  - local  stagnation  conditions  at  station  j 
r - atmospheric  conditions 

s - isentropic  conditions 

w - wall  conditions 

oo  - free  stream  conditions  for  flat -plate  flow 


INTRODUCTION 


A wind  tunnel  is  usually  considered  to  be  a device  in  which  models 
may  be  inserted  and  tested  under  controlled  conditions.  Wind  tunnels  have 
been  also  used,  however,  to  study  basic  phenomena  in  fluid  mechanics, 
such  as  boundary -layer  mechanics,  turbulence,  stability  of  the  laminar 
boundary  layer,  interaction  of  shock  waves  and  boundary  layers,  and,  in 
more  recent  years,  heat-transfer  phenomena  in  supersonic  flow.  Most 
supersonic  wind  tunnels  which  are  large  enough  to  insert  models  of  reason- 
able size  are  expensive  to  build  and  to  operate,  and  also  expensive  when 
used  to  measure:  heat-transfer  data  for  supersonic  flow.  Supersonic  flow 
established  under  controlled  conditions  in  a small  round  tube,  cne-half 
inch  in  diameter  can  be  used  to  measure  local  heat-transfer  coefficients. 
The  corresponding  apparatus,  which  is  described  herein,  may  be  considered 
to  be  an  example  of  an  inexpensive  supersonic  wind  tunnel,  in  that  it  is 
both  relatively  cheap  to  design  and  construct,  and  also,  relatively  low  in 
operational  costs. 


A recent  survey  (1)*  has  shown  the  great  *»eed  for  reliable  heat- 
Iransfer  data  for  supersonic  flow  in  view  of  the  small  amount  of  such  data 
published  during  the  years  from  1948  to  1953.  For  this  reason  a fairly 
detailed  account  of  the  results  obtained  in  an  investigation  of  heat-transfer 
coefficients  for  supersonic  flow  of  air  in  a round  tube  is  to  be  presented  in 
a series  cf  papers.  The  first  paper,  which  is  the  present  one,  consists  of 
a description  of  the  two  sets  of  apparatus  used,  the  original  data,  and  the 
calculated  results  for  a laminar  boundary  layer  based  on  a simple  one- 
dirnensionai  flow  model.  The  second  paper  will  present  the  calculated 
results  for  a laminar  boundary  layer  based  on  a two-dimensional  flow  model 
for  the  entrance  region  of  a tube.  The  third  paper  will  present  the  original 
data  and  the  calculated  results  for  a transitional  and  for  a turbulent  boundary 
layer  based  on  a simple  one -dimensional  flow  model. 

The  present  investigation  of  heat-transfer  coefficients  for  super- 
sonic flow  of  air  in  a round  tube  was  started  some  seven  years  ago  under 
the  sponsorship  of  the  Office  of  Naval  Research.  Preliminary  data  and 
results  were  obtained  several  years  ago  with  the  first  apparatus,  described 
herein  as  test  combination  C.  Similar  results  were  not  then  available  in 
the  literature  for  comparison.  The  possibility  of  systematic  errors  in 
the  data  and  the  lack  of  knowledge  of  the  type  of  boundary  layer  which 
existed  in  the  supersonic  flow  of  air  in  the  tube  were  two  main  reasons 
for  designing,  building,  and  testing  a new  apparatus  to  get  more  reliable 
heat-transfer  coefficients  for  supersonic  flow. 

♦Numbers  in  parentheses  refer  to  the  Bibliography  at  the  end  of  the  prper. 
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In  the  last  few  years  the  results  obtained  for  supersonic  flow  of 
air  in  a round  tube,  with  and  without  heat  transfer  to  the  air  stream, 
have  been  placed  on  a firm  foundation  as  indicated  by  the  following 

oin  Hpnrp  • 


1.  The  experimental  data  obtained  by  different  teams  of  students 
working  with  different  test  combinations,  for  both  adiabatic  and  diabatic 
flow,  have  shown  the  absence  of  significant  systematic  errors. 

2.  The  nature  of  the  boundary  layer  present  in  the  tube  flow  has 
been  conclusively  demonstrated  by  means  of  accurately  measured  velocity 
profiles  (2,3). 

3.  The  accurately  measured  velocity  profiles  for  supersonic  flow  in 
the  tube  have  also  demonstrated  the  basic  soundness  of  the  two-dimensional 
flow  model  used  to  interpret  the  data. 

4.  theoretical  solutions  (4)  have  been  obtained  for  the  system  of 
partial  differential  equations  of  energy,  momentum,  and  continuity  for  the 
case  of  supersonic  flow  of  air  in  a tube  in  the  entrance  region.  These  solu- 
tions describe  the  behaviour  of  supersonic  flow  in  a tube  with  a laminar 
boundary  layer  developing  from  the  entrance  plane  of  the  tube. 

METHODS  OF  TESTING 

The  experimental  program  described  here  has  consisted  of  two 
parts.  In  the  first  part,  a well-insulated  round  tube  was  used  to  measure 
values  of  the  local  adiabatic  wall  temperature  and  local  static  pressure 
of  a supersonic  stream  of  air.  The  results  of  these  tests,  with  a complete 
description  of  the  apparatus,  are  given  in  references  (5)  and  (6).  In  order 
to  minimize  the  amount  of  repetition,  the  reader  is  referred  to  these  two 
papers  for  many  details  not  given  in  this  paper.  In  the  second  part  of  the 
program,  two  different  test  combinations  were  used  in  which  steam  was 
condensed  outside  a round  brass  tube  in  order  to  measure  the  local  co- 
efficient of  heat  transfer  to  a supersonic  stream  of  air  flowing  inside  the 
tube. 


The  general  preparation  of  the  air  flow  and  the  method  of  drying 
the  air  stream  are  given  in  detail  in  reference  (5).  These  are  the  same 
fer  both  parts  of  the  program.  The  schematic  layout  of  the  entire  flow 
system  is  shown  in  Fig,  1.  The  purpose  of  the  major  changes  made  for  the 
tests  with  diabatic  flow  was  to  insure  that  almost  saturated  or  slightly 
superheated  steam  entered  the  outer  steam  chest  of  the  apparatus  and  to 
minimize  the  heat  losses  from  the  apparatus  to  the  surroundings. 
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A major  difference  between  tests  for  adiabatic  and  diabatie  flow 
was  the  length  of  time  required  to  bring  the  large  mass  of  heat-transfer 
apparatus  up  to  a fixed  temperature  corresponding  to  that  of  the  condensing 
steam.  This  transient  heating  process  was  carefully  observed  and  recorded 
in  the  diabatie.  tests  until  r steady-state  condition  was  obtained.  Heat- 
transfer  measurements  were  made  oniy  for  the  steady  state.  During  the 
large  time  interval  required  for  both  the  transient  process  and  the  meas- 
urements, the  conditions  in  the  steam  main  varied  only  slowly.  Control 
devices  were  found  necessary  to  superheat,  or  desuperheat  the  steam  only 
slightly  in  order  to  maintain  constant  the  state  of  the  steam  fed  to  the  steam 
chest  of  the  apparatus.  They  are  indicated  in  the  scnematic  layout  of 
Fig.  1.  The  measured  superheat  was  found  to  be  about  0.5°  F on  the 
average. 


EXPERIMENTAL  APPARATUS 

Some  salient  features  of  the  two  test  combinations  of  heai-transfer 
apparatus  are  described  in  Table  1.  The  length-diameter  ratio  in  Table  1 
is  based  on  the  distance  from  the  end  of  the  curved  contour  of  the  super- 
sonic nozzle  to  the  exit  plane  of  the  test  section. 

TABLE  1 


TEST  COMBINATIONS  OF  HEAT-TRANSFER  APPARATUS 


SYMBOL 

NOZZLE 

TEST  SECTION 

L/D 

M 

o 

C 

Brass 

Brass 

31.8 

3.0 

D 

Stainless  steel 

Brass 

49.8 

2.8 

The  nouzles  used  to  produce  supersonic  flow  in  test  combinations 
C and  D were  different  in  design  and  in  final  contour.  Details  of  design 
and  of  constri  ction  of  these  two  nozzles  are  shown  in  reference  (5).  For 
purposes  of  illustration,  the  rough  dimensions  and  shapes  of  these  two 
nozzles  are  shown  in  Fig.  2.  In  the  brass  nozzle  of  test  combination  C,  no 
means  of  measuring  the  temperature  gradient  along  the  axial  length  was 
available,  whereas  two  thermocouples  were  installed  in  the  stainless-steel 
nozzle  to  determine  the  axial  temperature  gradients  for  test  combination  D. 


Design  Considerations  for  Test  Sections 

The  first  xest  section,  used  in  test  combination  C was  designed  to 
produce  a small  thermal  resistance  for  radial  heat  flow  through  the  tube 
wall  in  comparison  with  the  thermal  resistance  being  measured,  namely, 
that  from  the  wall  to  the  air  stream  in  tne  tube.  Brass  was  therefore 


1 


i 


4 


selected  for  the  tube  material.  The  temperature  of  the  tube  was  maintained 
constant  over  its  entire  length  by  surrounding  it  with  condensing  steam. 

The  pvnprience  gained  after  operation  of  teBt  combination  C for 
about  two  years  resulted  in  the  following  additional  design  requii cments 
for  test  combination  D. 

1.  The  joint  between  the  nozzle  and  collar,  and  that  between  the 
collar  and  test  section  should  be  made  as  nearly  perfect  as  possible.  A 
misfit  in  excess  of  0.0002  in.  in  the  diameter  of  0.5  in.  is  undesirable. 

2.  The  axial  heat  flow  from  the  hot  test  section  to  the  cooler  nozzle 
should  be  reduced  to  a minimum.  This  objective  was  achieved  in  test  sec- 
tion D by  use  of  a nozzle  and  a collar  made  of  stainless  steel  which  has  a 
considerably  lower  thermal  conductivity  than  that  of  brass 

3.  The  no-load  or  blank-run  condensation  rate,  which  is  determined 
by  flow  of  steam  into  the  steam  chest,  for  zero  flow  of  air  in  the  test  section 
and  which  is  an  indication  of  the  total  amount  of  undesired  heat  flow  from 
the  steam  to  the  surroundings,  should  be  minimized,  especially  when  a 
laminar  boundary  layer  is  established  in  the  test  section. 

4.  The  state  of  the  steam  entering  the  steam  chest  must  be  care- 
fully controlled  to  be  slightly  superheated,  and  the  steam  flow  distribution 
into  and  within  the  steam  chest  must  be  as  nearly  uniform  as  practical. 


Test  Apparatus 

In  order  to  minimize  repetition,  the  general  common  features  of 
test  combinations  C and  "*  will  be  described  first.  The  dried  air  stream 
leaves  the  upstream  stagnation  tank,  passes  through  the  supersonic  nozzle 
into  the  half-inch  test  section  of  brass,  and  then  passes  out  through  the 
downstream  stagnation  rank  into  the  ejector  to  the  rtmosph.r 
drawn  from  the  supply  main  passes  through  a trap,  a pressure  regulator, 
a combination  of  a desuperheater  and  a superheater  with  electrical  heat 
input,  enters  the  steam  chest  through  the  inlet  pipe,  and  is  distributed 
around  the  test  section.  The  steam  is  prevented  from  flowing  directly  to 
the  outside  of  the  test  section  by  means  of  a brass  semicylindrical 
umbrella  which  covers  the  entire  length  of  the  tube.  This  umbrella  serves 
to  prevent  any  condensate  other  than  that  formed  on  the  tube  from  entering 
the  compartments  around  the  test  section.  It  serves  also  as  a radiation 
shield. 


The  outside  of  the  test  section  is  partly  surrounded  by  a serni- 
c>Andrical  brass  trough,  open  at  the  top.  This  trough  is  subdivided  into 
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compartments  by  means  of  thin  brar-^  ; irtiticns  set  perpendicular  to  the 
tube  axis  at  definite  intervals.  The  steam  in  contact  with  the  tube  wall 
in  a given  compartment  condenses  at  a rate  proportional  to  the  rate  of 
heat  tranr'er  to  the  air  flowing  inside  the  tube.  This  condensate  is  caught 
in  the  trough  and  drains  from  the  compartment  through  a short  piece  of 
neoprene  tubing  into  a steam -jacketed,  calibrated,  and  graduated  glass 
collecting  tube.  The  glass  collecting  tubes  are  shielded  from  the  atmos- 
phere by  highly  reflective  aluminum  foil. 

The  entire  test  section,  together  with  nozzle  and  stagnation  tanks, 
is  covered  with  a large  thickness  of  insulating  felt  and  then  with  aluminum 
foil  to  reduce  extraneous  heat  flows.  The  total  thermal  capacity  of  this 
insulation,  the  test  section,  and  the  surrounding  steam  chest  is  quite  large, 
so  that  a large  amount  of  steam  condenses  in  the  transient  warm-up  process 
before  the  steady-state  measurements  of  heat  transfer  are  made.  Drains 
for  removing  the  condensate  and  atmospheric  vents  to  speed  up  preliminary 
heating  were  provided  in  both  test  combinations. 

The  method  of  construction  of  the  test  section  was  almost  th  „ same 
for  the  two  test  combinations.  An  axial  hole  was  drilled  in  •?.  «c!id  brass 
rod  and  then  the  numerous  holes  for  pressure  taps  were  drilled  through  the 
wall  of  the  tube.  The  inside  surface  of  the  tube  was  prepared  by  successive 
polishing  operations  alternating  with  successive  cleanings  of  the  pressure - 
tap  holes,  until  a high  polish  was  present  on  the  inside  surface  with  no 
detectable  burrs  at  the  pressure  taps.  The  angular  position  of  the  pressure - 
tap  holes  was  rotated  from  the  entrance  of  the  tube  to  the  exit. 

The  steam  chesl  for  test  combination  C was  almost  square  in  cross 
section  and  was  made  from  welded  plates,  except  for  two  sides  c i iiO  chesi 
for  test  combination  D was  made  from  a circular  casting,  with  two  halves 
joined  by  bolts  at  a gasket,  in  order  to  be  able  to  increase  the  working 
pressure  of  the  steam  in  the  chest. 

The  measurements  made  for  both  test  combinations  C and  D included 
the  upstream  stagnation  temperature  and  pressure,  the  tube -wall  tempera- 
ture along  its  entire  length,  the  local  static  pressure  of  the  air  flow,  the 
local  rate  of  collection  of  condensate  from  each  compartment,  the  down- 
stream stagnation  temperature  and  pressure,  the  no-load  condensate  rate 
for  zero  air  flow,  and  the  pressure  and  temperature  of  the  steam  at  several 
locations  in  the  steam  jacket. 

Test  Combination  C 

The  details  of  test  combination  C are  shown  in  Fig.  3.  Photographs 
of  this  apparatus,  before  and  after  assembly,  are  shown  in  Figs.  4 and  5, 
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respectively.  The  compartments  of  the  condensate -collecting  trough  are 
clearly  visible  in  Fig.  4.  The  oreteeting  umbrella  is  clearly  visible  in 
Fig.  5,  which  also  shows  the  pressure  connections,  the  condensate  tubing, 
and  steam -jacketed  glass  collecting  tubes.  The  upstream  stagnation  tank 
is  to  the  left  in  Fig.  5.  The  important  dimensions  of  test  combination  C 
are  summarized  in  Appendix  A. 

When  test  combination  C wa3  first  assembled,  a hard  rubber  collar 
was  used  between  the  nozzle  and  the  test  section.  The  first  t welve  runs 
(K-l  through  K-12)  yielded  pressure  distributions  which  were  not  as  smooth 
as  those  obtained  for  the  measurements  with  adiabatic  flow.  The  nozzle  and 
collar  were  taken  apart  and  it  was  found  that  they  were  misaligned  by  about 
0.001  in.  A new  collar  was  prepared  and  the  joint  between  nozzle  arid  collar 
and  that  between  collar  and  test  section  were  carefully  polished  and  aligned. 
The  remaining  heat-transfer  runs  made  with  test  combination  C were  found 
to  yield  smooth  pressure  distributions,  indicating  the  achievement  of  smooth 
flow  at  the  entrance  of  the  test  section. 


Test  Combination  D 

The  details  of  test  combination  D are  shown  in  Fig.  6.  Photographs 
of  this  apparatus,  before  and  after  assembly,  are  shown  in  Figs.  7 and  8, 
respectively.  Fig.  7 shows  clearly  the  stainless-steel  supersonic  nozzle, 
the  collar,  the  test  section  with  the  surrounding  compartments,  the  large 
castings  for  the  steam  chest,  and,  finally,  the  downstream  s'  ^nation  tank. 
Fig.  8 shows  the  connections  for  the  pressure  leads,  for  thi  lermocouples 
in  the  test  section,  and  for  the  condensate  collecting  tubes. 

Test  combination  D differed  from  C mainly  in  the  use  of  a longer 
test  section.  The  length-diameter  ratio  for  D is  49.8  compared  to  31.8  for 
C.  This  longer  length  of  D was  used  purposely  to  check  the  prediction, 
based  on  the  results  of  adiabatic  flow  interpreted  by  means  of  the  two- 
dimensional  flow  model,  that  a much  longer  laminar  boundary  layer  could 
be  maintained,  even  with  heat  transfer  to  the  air  stream,  than  was  used  in 
test  combination  C.  The  results  given  here  for  test  combination  D confirm 
this  prediction  quite  well. 

Test  combination  D has  also  been  used  to  measure  the  velocity 
profiles  for  supersonic  flow  in  the  test  section,  with  and  without  heat  trans- 
fer to  the  air  stream.  The  preliminary  results  obtained  from  velocity 
profile  measurements,  discussed  in  references  (2)  and  (3),  confirm  that 
either  a laminar  boundary  layer  can  exist  for  the  entire  tube  length,  or 
that  transition  of  this  boundary  layer  to  a tui’bulent  one  may  occur  before 
tube  exit.  Thus  the  interpretation  of  ail  previous,and  present  data  obtained 
with  this  type  of  apparatus  is  placed  on  a secure  foundation. 
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EXPERIMENTAL  RESULTS 

Appendix  A contains  the  measurements  for  flow  with  heat  transfer 
obtained  with  test  combinations  C and  D,  together  with  the  calculated 
results  bas«_d  on  a simpl  e one -dimensional  flow  model.  In  order  to  keep 
the  size  of  this  paper  within  reason,  only  those  data  which  correspond 
mainly  to  a laminar  boundary  layer  existing  over  most  of  the  tube  length 
are  given  here;  the  remaining  data,  corresponding  to  a boundary  layer  in 
transition  and  to  a iurbuieni  boundary  layer  will  be  given  in  similar  de- 
tail in  a later  paper  of  this  series.  Appendix  B contains  the  analysis  for 
the  simple  one -dimensional  flow  model,  the  detailed  method  of  computation, 
and  a sample  calculation. 

The  number  of  runs  made  with  test  combinations  C and  D,  the 
ranges  of  Reynolds  numbers  covered  and  the  prevailing  type  of  boundary 
layer  present  in  the  tube,  are  given  in  Table  2. 


TABLE  2 


SUMMARY  OF  HEAT -TRANSFER  TESTS 


TEST 

TYPE  OF 

NO.  OF 

r"d 

Re_ 

COMBINATION 

BOUNDARY 

RUNS 

L 

LAYER 

MIN. 

MAX. 

MAX. 

C 

Laminar 

4 

52,000 

96,000 

2,270,000 

C 

f Turbulent 
l Transitional 

19 

43,000 

380,000 

8,560,000 

D 

Laminar 

13 

22,000 

92,000 

3,590,000 

D 

/ Turbulent 

L Transitional 

G 

20,000 

492,000 

12,000,000 

The  diameter  Reynolds  number  is  based  on  the  tube  diameter  and  on  the  mean 
stream  properties  measured  or  computed  at  the  first  station  — for  C the 
first  station  is  1.66  in,  from  the  exit  plane  of  the  nozzle  whereas  for  D it 
is  1.19  in.  from  the  exit  plane  of  the  nozzle.  Hence  the  inlet  diameter 
Reynolds  numbers  for  C do  not  correspond  exactly  to  those  for  D.  The 
length  Reynolds  number  is  based  on  the  distance  from  the  end  of  the  curved 
contour  of  each  nozzle  and  on  the  mean  stream  properties  at  that  distance. 
Hence  the  values  of  ReL  for  C and  D are  less  arbitrary  than  those  for 
ReD  and  more  nearly  comparable,  especially  at  a considerable  distance 
downstream  from  the  entrance  plane  of  the  tube. 
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Laminar  Boundary  Layer 

The  data  presented  here  for  a laminar  boundary  layer  in  a round 
tube  were  first  interpreted  on  the  basis  of  a simple  one -dimensional  flow 
model.  This  model  ignores  completely  the  growth  of  the  laminar  bounaary 
layer  in  the  tube  and  also  the  fact  that  this  type  of  supersonic  flow  is  mainly 
one  of  "entrance  flow"  with  insufficient  flow  length  to  produce  a "fully- 
developed"  flow.  Moreover,  previous  work  (5)  on  adiabatic  supersonic  flow 
in  a tube  has  shown  that  this  model  is  inadquate  for  calculation  of  friction 
coefficients  and  recovery  factors  for  a laminar  boundary  layer.  In  the 
present  paper  this  simple  model  is  still  used  since  it  permits  one  to  get 
a quick  reduction  of  the  original  data  to  useful  form,  but  the  computed 
quantities  arc  interpreted,  related,  and  compared  with  essentially  the  cor- 
rect phenomenological  picture  of  the  growth  of  a laminar  boundary  layer 
either  on  a flat  plate  or  in  a tube. 

For  flow  over  a flat  plate,  a laminar  boundary  layer  begins  to  grow 
from  the  leading  edge  until  it  undergoes  a transition  to  become  a turbulent 
boundary  layer.  The  transition  process  occurs  in  a finite  length  of  flow. 

The  turbulent  boundary  layer  continues  to  grow  in  thickness  in  the  direction 
of  flow.  This  simple  picture  of  plate  flow  can  be  used  with  precision  to 
develop  a two-  dimensional  flow  model  for  supersonic  flow  in  a tube  (6). 

The  details  of  this  two-dimensional  flow  model  for  tube  ilow  with  heat 
transfer  will  be  given  in  a later  publication.  The  results  obtained  with 
iiiis  more  exact  model  justify  the  method  of  comparison  of  tube  flow  with 
plate  flow  which  is  used  in  the  present  paper. 

The  results  computed  on  the  basis  of  the  one -dimensional  flow  model 
will  also  be  compared  with  the  theoretical  results  for  tube  flow  obtained 
recently  in  this  program  (7).  These  results  were  obtained  by  investigation 
of  the  basic  partial-differential  equations  of  energy,  momentum,  and 
continuity  for  a developing  laminar  boundary  layer  adjacent  to  an  isentropic 
core  in  the  central  portion  of  the  tube.  After  transformations,  these  equa- 
tions were  solved  with  the  aid  of  the  M.I.T.  Differential  Analyzer.  The 
solutions  were  obtained  on  the  basis  of  the  simplifying  assumption  of  con- 
stant fluid  viscosity  and  thermal  conductivity.  The  corresponding  cal- 
culated results  will  be  referred  to  on  the  charts  to  follow  as  "tube  flow, 
constant  <i  and  X." 

The  experimental  results  for  tube  flow  with  a laminar  boundary  layer, 
covering  the  range  of  inlet  diameter  Reynolds  number  from  about  20,000  to 
100,000,  are  presented  in  six  charts.  Figs.  9 to  14.  Each  chart  shows  the 
measured  values  of  the  modified  pressure  ratio,  the  wall  temperature  of 
the  brass  test  section,  the  heat  flux,  and  finally  the  computed  values  of  the 
Stanton  number . Each  of  these  quantities  is  plotted  against  the  length 
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Reynolds  number.  Ail  the  data  for  predominantly  laminar  boundary  layers 
which  were  obtained  with  test  combinations  C and  D are  given  on  these 
six  charts;  the  remaining  data  for  a predominantly  transitional  or  turbulent 
boundary  layer  are  given  in  a later  paper  and  contain  only  a small  amount 
of  information  on  the  laminar  boundary  layer.  Six  charts  were  found  to  be 
necessary  to  present  the  data  in  sufficient  detail  because  they  were  found 
to  be  quite  sensitive  to  small  changes  in  values  of  the  inlet  diameter 
Reynolds  number. 

The  ratio  of  measured  local  static  pressure  to  stagnation  pressure 
is  given  in  Figs.  3 to  14  in  terms  of  a dimensionless  modified  pressure 
ratio  mainly  to  place  results  from  test  combinations  C and  D and  runs 
made  with  different  stagnation  temperatures  on  a comparable  basis.  The 
local  wall  temperature  could  have  been  presented  in  terms  of  the  ratio  of 
local  wall  temperature  to  stagnation  temperature  but  the  actual  variations 
are  too  small  to  warrant  use  of  this  ratio.  The  local  heat  flux  is,  in 
reality,  an  average  value  of  heat  flux  over  the  short  length  of  one  of  the 
condensate  collecting  compartments. 

Fig.  9 presents  results  for  the  lowest  value  of  the  diameter  Reynolds 
number  atiaii.ed  in  these  heat-transfer  tests,  namely,  22,000.  Both  runs 
shown  in  Fig.  9 were  made  with  test  combination  D and  illustrate  the  degree 
of  reproducibility  of  the  data  for  static  pressure,  wall  temperature,  and 
heat  flux  along  the  length  of  the  test  section.  The  pressure  fluctuations 
shown  in  Fig.  9 are  not  uncommon  in  this  type  of  experiment.  The  wall 
temperature  is  remarkably  constant  along  the  length,  with  an  average  de- 
viation from  a mean  value  of  less  than  a fraction  of  one  degree.  The  heat 
flux  appears  to  be  smooth  in  the  entrance  portion  of  the  tube  which  cor- 
responds to  the  laminar  boundary  layer,  and  then  appears  to  fluctuate  about 
a mean  value  of  235  Btu/(hr  ft^). 

Fig.  9 compares  the  computed  Stanton  number  for  tube  flow  with  the 
values  predicted  for  plate  flow  and  for  tube  flow.  The  predicted  Stanton 
number  for  plate  flow  with  zero  pressure  gradient  is  taken  from  Van  Driest 
(8)  for  a laminar  boundary  layer  with  variable  viscosity  and  thermal  con- 
ductivity; the  comparison  is  made  for  a free-stream  Mach  number  ' 2,8 
and  a ratio  of  wall  temperature  to  free-stream  temperature  of  3.0,  The 
predicted  Stanton  number  for  tube  flow  is  taken  from  the  theoretical  solu  - 
tion given  by  Toong  (7)  for  a laminar  boundary  layer  with  constant  viscosity 
and  thermal  conductivity.  This  solution  shows  that  the  Stanton  number  de- 
pends both  on  the  length  Reynolds  number  and  the  diameter  Reynolds  number 
so  that  two  lines  are  shown  to  cover  the  range  of  diameter  Reynolds  number 
employed  in  the  tests.  For  the  entrance  region  of  the  tube,  where  a 
laminar  boundary  layer  is  forming,  the  experimental  tube  results  are  in 
good  agreement  with  the  theoretical  values  for  tube  flow.  Tht  experimental 
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data  lie  below  the  theoretical  values  for  plate  flow  probably  because  of 
the  pretence  of  a large  adverse  pressure  gradient  in  the  tube  flow. 

Fig.  10  shows  the  results  for  a diameter  Reynolds  number  of  32,000  for 
test  combination  D.  The  pressure  data  are  smoother  and  more  reproducible 
than  those  in  Fig.  9,  and  the  heat-flux  data  cover  a much  larger  range  of 
values  than  those  in  Fig.  9.  The  heat-flux  data  and  the  Stanton  numbers  hnth 
indicate  a sharp  rise  in  value  near  the  exit  of  the  test  section;  this  sharp 
rise  is  one  of  the  first  signs  noted  of  transition  from  a laminar  to  a turbulent 
boundary  layer  in  the  tube.  A f in  the  previous  case,  where  a laminar  boundary 
!ayer  is  present,  the  values  of  the  Stanton  number  compart  well  with  the 
theoretical  predictions  for  tube  flow  but  lie  below  those  for  plate  flow  with 
zero  pressure  gradient. 

Fig.  11  shows  the  results  for  3 runs  with  test  combination  D and 
one  run  for  test  combination  C.  The  sensitivity  of  the  data  to  small  changes 
in  diameter  Reynolds  number  is  evident  in  Fig.  11,  where  the  local  heat- 
flux  values  for  Runs  B-14  and  K-20  arc  considerably  larger  than  those  for 
Runs  B-2  and  B-4  even  though  the  diameter  Reynolds  number  is  increased 
only  from  45,000  to  51,000.  The  modified  pressure  ratio  for  Run  K-20  is 
slightly  smaller  than  those  for  the  three  runs  with  test  combination  D,  but 
the  same  slow  rise  in  pressure  with  a maximum  value  attained  near  the 
downstream  end  of  the  tube  is  evident  is  all  four  runs.  Careful  study  was 
made  of  the  conditions  at  this  value  of  the  diameter  Reynolds  number  be- 
cause preliminary  data  for  adiabatic  flow  had  indicated  that  for  this  value 
the  laminar  boundary  layer  occupied  most  of  the  length  of  the  test  section. 
Substantiation  of  these  preliminary  data  is  seen  in  runs  B-2  and  B-4  where 
a laminar  boundary  layer  was  found  to  exist  up  to  station  17,  with  a value 
of  length-diameter  ratio  of  45.3.  The  agreement  of  the  measured  Stanton 
numbers  for  tube  flow  with  the  theoretical  values  for  plate  flow  and  tube 
flow  is  excellent  on  the  basis  of  the  simple  one -dimensional  flow  model. 

Similar  excellent  agreement  was  found  in  reference  (5)  for  the  local  apparent 
er  coefficients  of  adiabatic  tube  flow  at  a diameter  Reynolds  number 
of  000.  Fig.  11  also  indicates  that  the  start  of  the  transition  from  a 
laminar  to  a turbulent  boundary  layer  occurs  at  a value  of  the  length 
Reynolds  number  of  900,000  for  Run  K-20  and  about  1,300,000  for  the  other 
three  runs.  These  values  are  in  good  agreement  with  those  for  transition 
on  a flat  plate. 

In  Fig.  9 the  value  of  the  Stanton  number  begins  to  level  off  at  a 
length  Reynolds  number  of  200,000,  and  fluctuations  about  a constant  value 
extend  to  a Reyno’ds  number  of  fct00,000.  In  Fig.  10  a similar  behavior  is 
noted  between  Reynolds  numbers  of  300,000  and  900,000,  and  in  Fig.  11 
between  500,000  and  1,000,000.  A possible  explanation  of  this  levelling  off 
and  fluctuation  of  values  of  the  Stanton  number  could  be  the  phenomenon  of 
separation  of  the  laminar  boundary  layer  under  the  influence  of  an  adverse 
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pressure  gradient  which  is  present  over  a long  length  of  tube  flow.  The 
fluctuations  could  be  partially  explained  in  terms  of  reattachment  of  the 
separated  laminar  boundary  layer.  Further  work  is  needed  to  test  this 
explanation. 

Fig.  12  presents  the  data  for  a diameter  Reynolds  number  of  about 

70.000  for  2 runs  with  test  combination  D and  one  run  for  test  combination 
C.  The  excellent  agreement  of  the  values  of  the  modified  pressure  ratio 
and  of  the  heat  flux  for  the  laminar  boundary  layer  for  the  two  completely 
different  test  combinations  indicates  the  absence  of  appreciable  systematic 
errors.  The  second  striking  feature  of  the  data  in  Fig.  12  is  the  steep  rise  in 
the  heat  flux,  and  likewise  in  the  Stanton  number,  when  the  laminar  boundary 
layer  undergoes  a transition  to  a turbulent  one.  The  older  and  shorter  test 
combination  C was  not  long  enougli  to  show  this  feature.  The  steep  rise  in 
heat  flux  ooserved  here  is  an  excellent  indicator  of  transition.  The  static 
pressure,  and  the  wall  temperature’ on  the  other  hand,  changed  but  slightly 
over  the  entire  eight -diameters  of  length  of  this  transition  process.  The 
values  of  the  measured  Stanton  numbers  in  Fig.  12  are  in  excellent  agree- 
ment with  those  predicted  for  piate  flow,  and  are  larger  than  those  for  tube 
flow  with  a laminar  boundary  layer.  The  better  agreement  with  predicted 
values  for  plate  flow  may  also  be  due  to  the  smaller  pressure  gradient 
present  in  the  data  in  Fig.  12  than  in  Figs.  9,  10,  and  11. 

Fig.  13  shows  the  data  for  a diameter  Reynolds  number  of  about 

80.000  for  Runs  B-l  and  B-9.  The  results  are  similar  to  those  noted  for 
Fig.  12,  especially  with  respect  to  the  steep  slope  of  the  curve  for  heat 
flux  during  transition.  It  should  be  noted  that  as  the  vaiue  of  the  diameter 
Reynolds  number  is  increased  imm  Fig.  9 on  to  Fig.  14,  the  value  of  the 
length  Reynolds  number  at  transition,  computed  on  the  basis  of  the  one- 
dimensional  flow  model  gradually  increased  from  about  900,000  to  about 
2,200,000.  A similar  shift  was  found  in  the  transition  for  adiabatic  tube 
flow.  The  agreement  between  measured  values  of  the  Stanton  number 

and  those  predicted  for  plate  flow  is  excellent  for  this  type  of  measurement 
in  view  of  published  data. 

Fig.  14  shows  the  data  for  Runs  B -7,  B-16,  K-i9,  and  K-23  for  the 
last  and  largest  value  of  the  diameter  Reynolds  num  er,  90,000,  which  was 
chosen  to  demonstrate  the  data  for  a predominantly  laminar  boundary  layer. 
The  agreement  between  the  data  for  the  two  different  test  combinations  is 
again  excellent.  The  steep  slopes  of  the  curves  for  heat  flux  and  for  Stanton 
number  are  again  emphasized  in  the  data  taken  with  the  longer  test  combina- 
tion D.  Figs.  13  and  14  indicate  smaller  values  of  the  adverse  pressure 
gradient  than  do  the  preceding  figures.  For  these  four  runs,  the  agreement 
between  measured  Stanton  numbers  and  plate -flow  values  for  zero  pressure 
gradient  is  excellent. 
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One -Dimensional  Flow  Model 

Previous  work  (5,  6)  has  shown  that  the  one-dimensional  flow  model. 
1-DFM,  used  for  the  present  heat -transfer  calculations  is  not  fully  adequate 
for  calculation  or  understanding  of  adiabatic  supersonic  flow  in  the  entrance 
region  of  a tube.  The  experimental  results  given  in  Figs.  9 to  14,  additional 
data  on  measured  velocity  profiles  for  supersonic  tube  flow  (2:,  3),  and, 
finally,  unpublished  data  on  measured  temperature  profiles  for  supersonic 
tube  flow  ail  show  that  the  simple  1-DFM  is  not  fully  adequate  to  explain 
the  heat-transfer  data  for  supersonic  tube  flow. 

In  place  of  using  tube-type  flow,  as  the  1-DFM  would  suggest,  the 
data  have  been  treated  here  as  though  the  laminar  boundary  layer  develops 
in  the  tube  in  the  same  way  as  a laminar  boundary  layer  develops  on  a flat 
plate.  A somewhat  more  complicated  flow  model,  the  two-dimensional  flow 
model,  2-DFM,  will  be  used  in  a subsequent  paper  to  give  support  to  this 
phenomenological  explanation  of  supersonic  tube  flow  in  the  entrance  region. 

In  addition,  this  2-DFM  applied  to  heat-transfer  data  for  tube  flow  will 
indicate  the  means  of  transforming  the  tube  data  for  quantitative  comparison  with 
heat-transfer  data  for  plate  flow.  Finally  this  2-DFM  will  be  used  to  snow  the 
means  of  combining  all  the  data  for  a laminar  boundary  layer  in  tub^  flow  so 
as  to  reduce  the  scattering  inherent  in  the  measurement  of  local  values  of  the 
heat  flux  over  short  increments  of  tube  length. 

CONCLUSIONS 

Reliable  data  on  heat-transfer  coefficients  to  air  flowing  at  super- 
sonic velocities  in  a round  tube  are  presented  here  for  the  case  of  a lamina’' 
boundary  layer.  The  agreement  found  between  measurements  made  with 
completely  independent  test  combinations  by  different  groups  of  students  pro- 
vides assurance  that  no  significant  systematic  errors  exist. 

The  data  for  tube  flow  are  computed  on  the  basis  of  a simple  one- 
dimensional  flow  model.  Since  this  model  is  not  fully  adequate  to  explain  or 
interpret  the  heat -transfer  data  for  supersonic  tube  flow,  those  data  are  inter- 
preted in  terms  of  a laminar  boundary  layer  which  begins  to  grow  at  tube 
entrance.  The  calculated  values  are  compared  with  the  theoretical  predictions 
for  a laminar  boundary  layer  developing  over  a flat  plate  with  zero  pressure 
gradient  and  with  the  theoretical  predictions  for  tube  flow  based  on  constant 
viscosity  and  constant  thermal  conductivity. 

The  measured  Stanton  numbers  agree  best  with  the  theoretical  values 
for  tube  flow  and  deviate  most  from  the  values  for  plate  flow  for  the  lowest 
value  of  the  diameter  Reynolds  number  attained  here,  namely,  22,000,  On  the 


13 


other  hand,  the  reverse  is  true  at  the  highest  value  of  the  diameter  Reynolds 
number,  namely,  90,000.  The  agreement  between  measured  Stanton  numbers 
and  flat-plate  values  improves  as  the  adverse  pressure  gradient  in  the  tube 
decreases,  that  is,  as  the  value  of  the  diameter  Reynolds  number  increases 
from  20,000  to  90,000;  the  agreement  at  the  highest  values  of  the  diameter 
Reynolds  number  is  excellent  for  this  type  of  measurement. 

The  transition  from  a laminar  to  a turbulent  boundary  layer  for 
supersonic  flow  in  a tube,  with  heat  transfer  is  most  easily  detectable  hy 
the  sudden  sharp  rise  of  the  value  of  heat  flux  at  the  transition  region,  even 
though  hardly  any  change  in  static  pressure  is  evident.  The  numerical  value 
of  the  length  Reynolds  number  at  inception  of  transition  is  in  excellent  agree- 
ment with  similar  published  data  for  a laminar  boundary  layer  on  a flat  plate. 

The  data  indicate  that  a fairly  great  length  of  supersonic  flow  can  be 
established  in  a tube  with  a laminar  boundary  layer  over  most  of  this  length. 
Hence  this  type  of  apparatus  represents  an  inexpensive  supersonic  wind 
tunnel  for  adiabatic  and  diabatic  flow. 

The  one -dimensional  flow  model  was  used  in  this  paper  to  reduce  the 
experimental  values  to  computed  quantities  in  a simple  and  quick  way.  The 
computed  quantities,  however,  are  interpreted,  related,  and  compared  with 
essentially  the  correct  phenomenological  picture  of  the  growth  of  a laminar 
boundary  layer  in  a tube.  The  same  original  data  will  be  treated  in  a later 
paper  by  a more  detailed  analysis  based  on  a two-dimensional  flow  model. 
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G.  M.  Ketchum,  W.  O.  Young,  B.  W.  Birmingham,  L,  J.  Ingolfsrud,  W.  S. 

Wu,  and  P.  Ashurkoff,  Mrs.  A.  R.  Walker  assisted  in  some  of  the  compu- 
tational work.  Professor  W.  K.  McAdams  contributed  to  the  design  of  the 
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by  the  Office  of  Naval  Research  of  the  United  States  Navy. 
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APPENDIX  A 

MEASUREMENTS  AND  CALCULATED  RESULTS 

This  appendix  contains  the  detailed  dimensions  of  ihe  two  test  com- 
binations, the  original  measurements  lor  diabatic  flow,  and  the  calculated 
results  based  on  the  simple  one -dimensional  flow  model.  The  values  of  the 
discharge  coefficients  for  the  supersonic  nozzles  and  the  properties  of  air 
are  given  in  reference  (5). 

Details  of  Nozzles  and  Test  Sections 

The  values  of  the  nozzle  throat  diameter  D*,  and  the  test  section 
diameter  D,  are  given  in  Table  3 for  the  two  test  combinations, 

TABLE  3 


DIAMETERS  OF  NOZZLES  AND  TEST  SECTIONS 


TEST  COMBINATION 

NOZZLE 

TEST  SECTION 

D*  (in.) 

D (in.) 

C 

Brass 

Brass 

0,2416 

0.502 

D 

Stainless 

Steel 

Brass 

0. 2685 

0.5018 

The  length-diameter  ratios  of  the  two  test  sections  are  given  in 
Table  4.  The  length  is  the  distance  from  the  end  of  the  curved  contour 
in  the  supersonic  nozzle  to  the  location  of  the  pressure  tap  of  the  station 
in  question. 

The  value  of  the  local  heat-transfer  area  in  any  collecting  compart- 
ment is  found  by  multiplying  its  length  by  the  inner  circumference  of  the 
test  section.  For  tesi  combination  C-  the  length  of  each  collecting  compart- 
ment for  stations  1 through  14  is  1 in.,  but  for  the  first  collecting  compart- 
ment is  0.656  in.  (see  Fig.  3).  In  test  combination  D the  length  of  the 
compartments  for  stations  3 and  18  is  0.75  in.,  for  stations  4 through  9 
is  1.00  in.,  and  for  stations  10  through  17  is  2. GO  in.  (see  Fig.  6). 


Data  and  Calculated  Results 

The  original  data  and  calculated  results  for  seventeen  runs  corres- 
ponding to  flow  with  a laminar  boundary  layer  are  summarized  in  Table  5. 
Thirteen  runs  were  made  with  combination  D and  four  with  combination  C. 
The  data  include  measured  values  of  the  stagnation  temperature,  stagnation 


TABLE 


16 


f 


.O, 

< 

f- 


I 

I 

j 


Q 

Q 

Z 

< 

U 


c o 
Z 

o 

C-4 

<3 

Z 

w 

* £ 
o 
o 


w 

05 

co 

co 

W 

« 

0, 


H 

<3 

8 

J 


M 

O 

CQ 

e 

o 

• H 

k.  o 
nJ  a> 

3 w 
o +> 

U cn 
fU 


c 

Q5 

s 

t 

rt 

a 

E 

o 

U 

% 

Sh 


H Ck 


XJ 
C 
O 
u 
a> 
c o 


x 

"E  x x 

2 
H 


+->  r* 

ri  " 


u x x g s 5 

3 £ 3 > t»  q 

O -h  ,3  ni  ■ *h  • —3 

Ok  fa  co  co  W Z 


X 

■g 

X 

c 

J2 

<u 

at 

cu 

a> 

-*-» 

a> 

r~*. 

i 

E 

2 

l. 

3 

O 

a; 

-4-» 

<4-1 

• •H 

H 

H 

fa 

Ok 

a 

j 


CO 

«o  o»  t-  *r  r-  © 


oo 

CD 

CO 

[> 

t- 

CO 

CD 

lO 

<* 

CM 

*— t 

<X 

CD 

m 

»— 4 

CD 

CO 

■'t* 

CO 

CO 

CO 

co 

CO 

c 

• 

o 

iH 

CM 

CD 

CJO 

o 

03 

t- 

m 

03 

CO 

r- 

T— -1 

u5 

►—1 

r-i 

. i 

^Vl 

CM 

CM 

CO 

CO 

r Q, 
<1 


M 

W 

H 

Ok 

O 

co 

6 

HH 

CO 

55 

W 

r—H 

Q 


U 


W 

05 

D 

co 

CO 

W 

05 

CL 

Ok 

O 

o 

HH 

H 

<3 

U 

o 

J 


c 

a) 

E 

rt  5 

a 

s 

o 

U 

T3 

C 

O 

O 

a> 

CO 


T3 

U 

• H 

x 

H 


x: 

+* 

c 

0) 

> 


+■<  ^ x* 


X 

W) 


5 5 

Ok  co  co  Ok  fa 


X 

-+** 

c 

U> 

> 

(U 


C 

a> 

a> 

-4-» 
• H 

X 


M 4h 


f W E"1  E-* 


a> 

03 

L 

3 

O 

Ok 


fa 


i 


P 

/T\ 

o> 

T- < 

o> 

*0 

CD 

•o 

t- 

CD 

in 

Tt< 

CO 

CM 

*— 4 

o 

o 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

m 

r- 

05 

CO 

m 

C— 

05 

*H 

CO 

in 

r* 

05 

t— ^ 

*— < 

CM 

CM 

CM 

CM 

CM 

Z 

o 

I— I 

H 

<3 

H 

co 


d 

z 


N n ^ K!  ffl 


coroo^cvjco-^mcDr-co 


I 

LamK, 


48.08  Sixteenth 


17 


I 

i 


pressure,  local  wall  temperature,  local  static  pressure,  and  local  gross 
and  local  no  load  heat- transfer  rates.  The  calculated  results  are  based 
on  the  simple  one -dimensional  flow  model. 

Due  to  experimental  difficulties  during  certain  runs,  it  was  occasion- 
ally necessary  to  evaluate  certain  quantities,  such  as  wall  temperature, 
wall  pressure,  local  heat-transfer  rate,  by  linear  interpolation  or  extra- 
polation of  the  remaining  data.  Such  quantities  are  placed  in  parentheses 
in  Table  5. 

The  no-load  heat-iransfer  rates  were  not  determined  after  every  run. 
For  those  runs  where  such  data  were  available,  they  were  averaged  and  the 
average  values  were  used  also  for  the  runs  for  which  such  measurements 
were  omitted. 


Inspection  of  Table  5 shows  that  heat-transfer  parameters  are 
omitted  for  test  combination  D at  stations  3 and  18.  Actually,  the  net  heat- 
transfer  rate  measured  at  station  3 involves  the  additional  heat  transfer  to 
the  air  stream  which  occurs  in  the  nozzle,  collar,  and  a 0.50-inch  length 
of  test  section  upstream  of  station  3.  Similarly  at  station  18  the  net  heat- 
transfer  rate  involves  additional  heat  transfer  to  the  air  stream  in  a 0. 50- 
inch  length  of  test  section  downstream  of  station  18  and  through  the  down- 
stream boss  and  end  plate  to  the  air  in  the  downstream  stagnation  tank. 

The  data  at  these  two  stations  were  omitted  since  the/  cannot  be  corrected 
for  these  additional  heat  transfers.  The  analytical  results  of  Chapman  and 
Eubesin  (9),  together  with  measured  temperature  distributions  in  the  super- 
sonic nozzle,  were  used  to  check  the  measured  heat-transfer  rates  at  station 
3.  The  agreement  was  within  the  experimental  eri’or  of  the  measured  heat- 
transfer  rates. 


In  test  combination  C,  the  condensate  data  at  the  first  compartment, 
q0,  were  not  recorded.  The  amount  of  heat  transfer  to  the  air  stream  up- 
stream of  station  1 was  estimated  by  linear  extrapolation  of  the  data  from 
stations  i and  2. 
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APPENDIX  B 

ANALYSIS  AND  SAMPLE  CALCULATION 


Analysis 

Consider  the  steady  flow  of  air  from  a stagnation  state  through  a 
supersonic  nozzle  and  then  through  a round  tube  of  constant  cross-sectional 
area  in  which  heat  is  oorieri  to  the  air  stream.  The  following  assumptions 
are  made  in  the  analysis: 

1.  The  flow  is  one -dimensional,  i.e.  all  fluid  properties  are  uniform 
at  any  cross  section, 

2.  Air  is  a perfect  gas  with  a constant  value  of  the  ratio  of  specific 
heats  (k  = 1.40)  over  the  range  of  temperature  under  consideration. 

2,  Heat  added  to  the  air  stream  is  measured  by  the  amount  of 
condensate  collected  in  the  various  compartments. 


The  following  relations  hold  at  each  section  in  the  tube: 


Continuity: 

w = pVA 

(1) 

Equation  of  State: 

p = pRT 

(2) 

Energy: 

c_T  . = c T + Y“/2g 

P o..  p 

(3) 

Definition: 

M2  = V2/  gkKT 

(4) 

The  discharge  coefficient  of  the  supersonic  nozzle  is  defined  by 

c = (w/A*)/(w/ A*)  (5) 

w s 

For  isentropic  flow  to  the  nozzle  throat 


G*  = (w/A*) 
s s 


/ gk  k t 1 p* 

J 2 ,rT~’ 

' T Ol 


P*/P 


Ol. 


k l i 


(6) 


(7) 


For  a control  volume  which  encloses  the  fluid  between  the  upstream 
stagnation  slate  and  the  section  at  the  center  of  the  measuring  compartment 
of  station  j,  the  energy  equation  becorr.es  for  test  combination  C: 
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b 


i \ 

• l 
f 
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l 
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J 


-j-i 


q . / a 

J 


n = 0 


q = C (T  - T .). 
n p oi  oi 


(8a) 


For  test  combination  D the  energy  equation  is: 


n = i - 1 

q./2  + q = c:  (T  . - T ). 

J * n P oj  01 

1 1 - . ► 


(8b) 


Combining  Equations  (1)  through  (7), 


1 p A 

c n . A* 
w oi 


1 


K (-*_> 

IT.  k + 1 


_k_  + 1_ 

2 ( k 1) 


M 


f (M,k). 


. .2  1 
1VI 


(9) 


The  right-hand  side  of  F.quation  »9)  is  a unique  function  of  the  Mach  number 

if  k is  constant.  This  function  is  tabulated  under  the  heading  pA/p0A*  in 

Table  3G  of  reference  (10).  The  symbol  p ^ of  reference  (10)  identical 

with  n used  here, 
oi 

All  the  quantities  on  the  left-hand  side  of  Equation  (9)  are  measured 

except  cw  and  TQ-.  Equation  (8)  permits  the  calculation  of  TQ-  from 

measured  data.  The  plot  of  r versus  (Re*  ) , given  in  reference  (5), 

, vv  Li  s 

where 


(Fte*)  = G*  D*/p*  (10) 

D s s 

was  used.  Note  that  G|  can  be  found  from  Equations  (6)  and  (8)  from  measured 
data  and  that  p*  can  be  found  from  the  temperature. 

T*  = 0.83333T  = 0.83333<t  + 459.69).  (11) 

oi  oi 

\AAth  the  Mach  number  known,  the  mean-stream  temperature  can  be 
found  from 

k - 1 '* 

T /’v  _ 1/(1  TvT)  (12) 

m ' oj  ‘ 2 

which  is  tabulated  in  Table  30  of  reference  (10)  as  T/T  . The  following 
equations  were  used  to  complete  the  calculations: 


V 
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G 

= c G*  A*/ A 

(13) 

w s 

He.. 

= G D/m 

(1  4) 

JJ 

m 

Re, 

= GL in  = Re  (L/D) 

(15) 

Li 

m D 

h 

S,i/A'<T  -T  ) 

(16) 

w aw 

Nu 

; hD/A 

(17) 

D 

m 

Nu, 

= hL/A  = Nu(L/D) 

(18) 

m D 

St 

= h/c  G 

(19) 

P 

in  order  to  calculate  the  local  heat-transfer  coefficient  from  Equation 

(16),  It  is  necessary  to  calculate  the  adiabatic- wail  temperature  T . Fcr 

aw 

this  calculation  T was  found  lrom 

cl  W 


(20) 


In  order  to  calculate  the  local  heat -transfer  coefficients  from  Equa- 
tion (16),  the  adiabatic -wall  temperature  can  be  calculated  from  the  recovery 
factors  for  adiabatic  tube  flow  given  in  reference  (5).  Since  the  original 
data  on  which  these  recovery  factors  are  based  were  also  available  in  the 
form  of  the  i atio  of  adiabatic  -wall  temperature  to  stagnation  temperature, 
these  original  data  for  a laminar  boundary  layer  were  plotted  and  then- 
averaged.  An  average  value  of  (T„w/T0j)  of  0.9'10  was  used  in  this  paper 
In  order  to  reduce  the  time  required  for  these  computations. 


Sample  Calculation 

A sample  calculation  of  Run  No.  B-7  is  given  in  Table  6.  The  calcu- 
lations and  equation  numbers  given  in  the  headings  refer  to  the  preceding 
analysis.  The  values  given  in  Tables  5 and  6 are  based  on  calculations  made 
with  five  ^ r six  significant  figures  throughout. 
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SCHEMATIC  LAYOUT  OF  TEST  COMBINATIONS  C AND 
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fig.  3 detail:;  of  test  combination  c 


TEST  COMBINATION  C,  DEE  ORE  ASSEMBLY 
FIG.  4 


TEST  COMBINATION  C,  AFTER  ASSEMBLE 


FIG.  > 


Mifa^^{?li'fi^3}!SB;fc w--<  ffm  r **r  qfifl? 
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FIG.  6 DETAILS  OF  TEST  COMBINATION 
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TEST  COMBINATION  D,  BEFORE  ASSEMBLY 

FIG.  7 


TEST  COMBINATION  D,  AFTER  ASSEMBLY 
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